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„_jl We report results of a search for Weakly Interacting Massive Particles (WIMPs) with the Si 

detectors of the CDMS II experiment. This report describes a blind analysis of the first data taken 
with CDMS IPs full complement of detectors in 2006-2007. Results from this exposure using the 
Ge detectors have already been presented. We observed no candidate WIMP-scattering events in 
• • an exposure of 55.9 kg-days before analysis cuts. These data set an upper limit of 1.7 x 10 -41 cm 2 

on the WIMP-nucleon spin-independent cross section of a 10 GeV/c 2 WIMP; this limit improves 
to 8.3 x 10 -42 cm 2 in combination with previous Si data from this installation. These data exclude 
parameter space for spin-independent WIMP-nucleon elastic scattering that is relevant to recent 
searches for low-mass WIMPs. 

PACS numbers: 14.80.Ly, 95.35. +d, 95.30.Cq, 95.30.-k, 85.25.Oj, 29.40.Wk 



There is now overwhelming evidence that the bulk of 
the matter in our universe is in some nonluminous, non- 
baryonic form pQ. Though there is broad consensus on 
the amount of this dark matter present in the cosmos, 
its composition has thus far eluded laboratory investi- 
gations. Weakly Interacting Massive Particles (WIMPs) 
[2] - particles with masses between tens of GeV/c 2 and a 
few TeV/c 2 and interaction strengths characteristic of the 
weak force - form a leading class of candidates for this 



dark matter. Particles of this type would be produced 
thermally in the early universe in roughly the correct 
amount, and are predicted by many theoretical exten- 
sions to the Standard Model of particle physics [Tj, [3] . If 
WIMPs do constitute the dark matter in our galaxy, they 
may be detectable through their elastic scattering off of 
nuclei in terrestrial particle detectors @]. Numerous ex- 
perimental groups have sought to detect such scattering 
events using a wide variety of technologies [5]. 



2 



The Cryogenic Dark Matter Search (CDMS) collabora- 
tion seeks to identify nuclear recoils from WIMP interac- 
tions using semiconductor detectors operated at very low 
temperatures (—40 mK). These detectors use a simulta- 
neous measurement of ionization and out-of-equilibrium 
phonons to identify such events among a far more nu- 
merous background of electron recoils. From 2003-2008 
the collaboration operated CDMS II, an array of Ge and 
Si detectors located at the Soudan Underground Labora- 
tory. Previous results from the CDMS II installation j6j- 
[S] have set stringent upper limits on the WIMP-nucleon 
scattering cross section and constrained some non-WIMP 
dark matter candidates [TU] , 

This work presents results from a search for WIMP in- 
teractions in the CDMS II Si detectors during the first 
run of the experiment with its full complement of detec- 
tors. The lower atomic mass of Si generally makes it a less 
sensitive target for spin-independent (scalar) WIMP in- 
teractions, due to the coherent enhancement of the scat- 
tering cross section for heavy nuclei. The lower atomic 
mass of Si is advantageous in searches for WIMPs of rel- 
atively low mass, however, due to more favorable scat- 
tering kinematics. A WIMP of mass < 40 GeV/c 2 will 
impart more recoil energy to a Si atom than a Ge atom 
on average, so a WIMP of sufficiently low mass (M < 10 
GeV/c 2 for CDMS II) will generate more detectable re- 
coils in a Si detector at fixed energy threshold. New 
particles at such masses are generally disfavored in fits 
of supersymmetry models to precision electroweak data 
(e.g. [H]), but viable models in this regime do exist (e.g. 
[12]). Renewed interest in this mass range has been moti- 
vated by recent results from the DAMA/LIBRA (T3], Co- 
GeNT Q3], and CRESST [T5] experiments, which have 
been interpreted as possible evidence of WIMP scatter- 
ing. CDMS has previously explored similar parameter 
space using dedicated low-threshold analyses of data from 
its shallow and deep runs |16) . 

In its final configuration, the CDMS II array consisted 
of 30 Z-sensitive ionization and phonon (ZIP) detectors: 
19 Ge (-239 g each) and 11 Si (-106 g each), for a total of 
-4.6 kg of Ge and -1.2 kg of Si. Each CDMS detector is 
a semiconductor disk, 7.6 cm in diameter and 1 cm thick, 
instrumented to detect the phonons and ionization gener- 
ated by particle impacts within the crystal. One flat face 
of each detector is instrumented with four readout chan- 
nels composed of superconducting transition-edge sensors 
(TESs) to detect out-of-equilibrium phonons. The op- 
posite flat face is divided into two concentric ionization 
electrodes: an inner (primary) electrode covering —85% 
of the detector surface and an outer guard ring. The lat- 
ter defines a fiducial volume within each ZIP by identify- 
ing interactions near the detector rim, which may suffer 
from reduced ionization collection. We discriminate nu- 
clear recoils from background electron recoils using the 
ratio of ionization to phonon recoil energy ("ionization 
yield"), with a misidentification rate for electron recoils 
of less than 1 in 10 4 . Electron recoils that occur within 
—10 /xm of a detector surface can be misclassified as nu- 



clear recoils due to reduced ionization collection. Such 
surface events are identified by the faster arrival of their 
phonon signals, giving an overall misidentification rate 
for electron recoils of less than 1 in 10 6 . 

This detector array was housed within a low- 
radioactivity cryogenic installation [SJ [T7] at the Soudan 
Underground Laboratory, Minnesota, U.S.A. The rock 
overburden above the Soudan facility (2090 meters water 
equivalent) reduces the flux of cosmogenic muons inci- 
dent upon the detector installation by a factor of — 10 5 , 
thus greatly reducing the background neutron flux. An 
outer hermetic layer of plastic scintillator identifies re- 
maining cosmogenic muons entering the passive Pb and 
polyethylene shielding surrounding the detector volume. 

We consider data taken with the Si detectors during 
the first two run periods of the full CDMS II detector in- 
stallation, acquired between October 2006 and July 2007. 
The Ge results from this data set were described in a 
previous publication [8j, which was released before the 
Si analysis was complete. Of the 11 Si detectors, five 
were excluded from the WIMP-search analysis: two due 
to wiring failures that led to incomplete collection of the 
ionization signal, one due to unstable response on one 
of its four phonon channels, and two due to inadequate 
rejection of surface events in the analysis chosen for this 
exposure. Periods of poor performance, as identified by 
a series of Kolmogorov-Smirnov tests, were also excluded 
from analysis. After all such exclusions, these data rep- 
resent 55.9 kg-days of exposure with the remaining six Si 
detectors before selection of WIMP candidates. 

The response of these detectors to electron recoils was 
calibrated using events from extensive exposures to ra- 
dioactive 133 Ba sources in situ at Soudan. Electron re- 
coils from these sources were used to empirically charac- 
terize and correct for the dependence of phonon pulse 
shape on event position and energy. As in the anal- 
ysis of the Ge detectors, events at large detector radii 
were excluded due to degraded performance of this cor- 
rection technique. Because the Si detectors generally do 
not show a clear 356-keV spectral line from the 133 Ba 
source, their ionization and phonon energy scales were 
calibrated using 356-keV events that share their energy 
with a neighboring detector. 

The detectors' response to nuclear recoils was char- 
acterized using neutron-scattering events from in situ 
calibration with a 252 Cf source. The resulting nuclear 
recoil population was used to tune the various WIMP- 
selection criteria of this analysis, notably those for ion- 
ization yield and phonon timing. We have verified the 
calibration of the nuclear recoil energy scale by compar- 
isons to Monte Carlo simulations of the 252 Cf exposures, 
an analysis which will be described in a separate pub- 
lication [IS]. Such comparisons are particularly robust 
for the Si detectors, due to a fortuitous resonant feature 
in the Si- neutron elastic scattering cross section. This 
study indicates that our reconstructed energy may be 
10% lower than the true recoil energy in the relevant 
energy range. This would slightly weaken our quoted 
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FIG. 1. Nuclear recoil acceptance as a function of recoil en- 
ergy after successive application of each WIMP-selection cri- 
terion shown. The bold solid curve shows the overall efficiency 
of this analysis. The abrupt drops in acceptance at low re- 
coil energies reflect the elevated energy thresholds chosen for 
some detectors. 



results, as described below. 

Candidate WIMP-scattering events were identified by 
a series of selection criteria. These criteria were defined 
in parallel with those described in [5] for the Ge detectors 
using the same techniques. As with the Ge detectors, all 
WIMP-selection criteria were defined blindly using cal- 
ibration and masked WIMP-search data; for the latter, 
events in and near the WIMP-candidate region were au- 
tomatically masked from the data set during analysis and 
thus had no impact on the definition of the selection crite- 
ria. A WIMP candidate was required to have phonon and 
ionization signals inconsistent with noise alone, to exhibit 
no coincident energy in the scintillating veto shield or in 
any of the other 29 ZIP detectors, and not to be coinci- 
dent with beam spills of the NuMI neutrino beam [!§] . 
We further demanded that any candidate event occur 
within the detector's fiducial volume and have ionization 
yield and phonon pulse timing consistent with a nuclear 
recoil. The recoil energy of each candidate event must 
also lie below 100 keV and above a detector-dependent 
threshold ranging from 7 to 15 keV, also chosen blindly 
based upon calibration data. Fig. [I] shows the estimated 
fraction of WIMP-scatter events that would be accepted 
by these signal criteria. Signal acceptance was measured 
using nuclear recoils from 252 Cf calibration. Monte Carlo 
simulations indicate that multiple-scattered neutrons in 
calibration data reduce the measured efficiency of the 
fiducial volume selection by ^5.5% with respect to the 
true value for single-scatter nuclear recoils, so we have 
scaled its efficiency upward by this amount. Signal accep- 
tance is ^40% at most recoil energies, somewhat higher 
than that of the Ge analysis. After applying all selec- 
tion criteria, the exposure of this analysis is equivalent 
to 19.7 kg-days over a recoil energy range of 7-100 keV 
for a WIMP of mass 60 GeV/c 2 . 

Neutrons from cosmogenic or radioactive processes can 
produce nuclear recoils that are indistinguishable from 
those from an incident WIMP. Simulations of the rates 
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FIG. 2. Ionization yield versus recoil energy in all detectors 
included in this analysis for events passing all signal criteria 
except (top) and including (bottom) the phonon timing crite- 
rion. The curved lines indicate the signal region (±2cr from 
mean nuclear recoil yield) between 7 and 100 keV recoil en- 
ergies. Electron recoils in the detector bulk have yield near 
unity, above the vertical scale limits. 



of these processes using GEANT4 and FLUKA lead us 
to expect < 0.1 false candidate events in the Si detectors 
from neutrons in this exposure. 

A greater source of background is the misidentifica- 
tion of surface electron recoils, which may suffer from 
reduced ionization yield. As in the Ge analysis, we de- 
veloped a Bayesian estimate of the rate of misidentified 
surface events based upon the observed performance of 
the phonon timing cut for events near the WIMP-search 
signal region |20) . For the Si analysis we based our model 
only upon multiple-scatter events within the ionization 
yield acceptance region, since other event samples in- 
corporated into the Ge analysis were found to be less 
reliable predictors for Si. This model is not applicable 
to detectors at the top and bottom of their respective 
stacks, since it is impossible to identify multiple-scatter 
events on the outside face of such detectors. We thus 
decided to exclude detectors in these positions from this 
blind analysis. The final model predicts an average of 
l.l^Q g (stat.) ± 0.1(syst.) misidentified surface events in 
the six Si detectors during this exposure. 

After all WIMP-selection criteria were defined and the 
background estimate finalized, the signal regions of the 
Si detectors were unmasked on December 3, 2008. No 
candidate WIMP-scattering events were observed. Fig. [2] 
illustrates the distribution of events in and near the sig- 
nal region of the WIMP-search data set before (top) and 
after ( bottom) application of the phonon timing criterion. 
Fig. [3] shows an alternate view of these events, expressed 
in "normalized" versions of yield and timing that are 
transformed so that the WIMP acceptance regions of all 
detectors coincide. 

This null result constrains the available parameter 
space of WIMP dark matter models. We compute up- 
per limits on the WIMP-nucleon scattering cross section 
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FIG. 3. Normalized ionization yield (standard deviations 
from the nuclear recoil band centroid) versus normalized 
phonon timing parameter (/js from the timing criterion) for 
events in all detectors from the WIMP-search data set pass- 
ing all other selection criteria. The black box indicates the 
WIMP candidate selection region. Also plotted are nuclear 
recoils from 252 Cf calibration data (light, green dots). 
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FIG. 4. Comparison of 90% C.L. upper limits from these data 
(solid) with those from CDMS II Ge 0[]2] (dash, +), EDEL- 
WEISS [21] (x), XENON10 [22] (S2-only analysis, 0), and 
XENON100 [23] (□). The filled regions identify regions of in- 
terest associated with data from DAM A/LIBRA [HUM] (dark 
grey, 99.7% C.L.), CoGeNT [25] as interpreted by Kelso et al. 
|26| (magenta, 90% C.L., including the effect of a residual 
surface event contamination), and CRESST II [15] (yellow, 
95.4% C.L.). 

using Yellin's optimum interval method |27| : this is equiv- 
alent to a Poisson upper limit in the present zero-event 



case, but generally results in a stronger limit when events 
are observed. We work within the "standard" halo model 
described in [55], assuming a galactic escape velocity of 
544 km/s [55] • Fig. [| shows upper limits on the WIMP- 
nucleon spin-independent scattering cross section at the 
90% confidence level from CDMS II data and a selection 
of other recent results. The present data set an upper 
limit of 1.66 x 1CT 41 (1.86 x 1CT 42 ) cm 2 for a WIMP 
of mass 10 (60) GeV/c 2 . The effect of a -10% increase 
in our nuclear recoil energy scale is well approximated 
below 20 GeV/c 2 by shifting the limit curve parallel to 
the mass axis by —7%. Since unblinding these data, re- 
cent results from CDMS II US], EDELWEISS [2T], 
XENON100 [23], and a novel low-threshold analysis of 
data from XENON10 [55J also disfavor this parameter 
space. 

Fig. [4] also compares these results to three recent re- 
sults that have been interpreted as evidence for WIMP 
interactions. The CoGeNT experiment has observed an 
excess of events in their Ge crystal above expected back- 
ground [14] and an annual modulation of their low-energy 
event rate [25] [26], similar to what might be expected 
from interactions of a low-mass WIMP. The CRESST II 
experiment has also observed an excess of events above 
their background model [15]. This null result disfa- 
vors portions of the best-fit regions suggested by the 
authors in both cases, as well as an interpretation of 
the DAMA/LIBRA annual modulation signal in terms 
of spin-independent scattering [24] . 
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